Universal method for the isolation of microvessels from frozen brain tissue: A proof-of-concept multiomic investigation of the neurovasculature.

The neurovascular unit, comprised of vascular cell types that collectively regulate cerebral blood flow to meet the needs of coupled neurons, is paramount for the proper function of the central nervous system. The neurovascular unit gatekeeps blood-brain barrier properties, which experiences impairment in several central nervous system diseases associated with neuroinflammation and contributes to pathogenesis. To better understand function and dysfunction at the neurovascular unit and how it may confer inflammatory processes within the brain, isolation and characterization of the neurovascular unit is needed. Here, we describe a singular, standardized protocol to enrich and isolate microvessels from archived snap-frozen human and frozen mouse cerebral cortex using mechanical homogenization and centrifugation-separation that preserves the structural integrity and multicellular composition of microvessel fragments. For the first time, microvessels are isolated from postmortem ventromedial prefrontal cortex tissue and are comprehensively investigated as a structural unit using both RNA sequencing and Liquid Chromatography with tandem mass spectrometry (LC-MS/MS). Both the transcriptome and proteome are obtained and compared, demonstrating that the isolated brain microvessel is a robust model for the NVU and can be used to generate highly informative datasets in both physiological and disease contexts.

SNORD90 induces glutamatergic signaling following treatment with monoaminergic antidepressants.

Pharmacotherapies for the treatment of major depressive disorder were serendipitously discovered almost seven decades ago. From this discovery, scientists pinpointed the monoaminergic system as the primary target associated with symptom alleviation. As a result, most antidepressants have been engineered to act on the monoaminergic system more selectively, primarily on serotonin, in an effort to increase treatment response and reduce unfavorable side effects. However, slow and inconsistent clinical responses continue to be observed with these available treatments. Recent findings point to the glutamatergic system as a target for rapid acting antidepressants. Investigating different cohorts of depressed individuals treated with serotonergic and other monoaminergic antidepressants, we found that the expression of a small nucleolar RNA, SNORD90, was elevated following treatment response. When we increased Snord90 levels in the mouse anterior cingulate cortex (ACC), a brain region regulating mood responses, we observed antidepressive-like behaviors. We identified neuregulin 3 (NRG3) as one of the targets of SNORD90, which we show is regulated through the accumulation of N6-methyladenosine modifications leading to YTHDF2-mediated RNA decay. We further demonstrate that a decrease in NRG3 expression resulted in increased glutamatergic release in the mouse ACC. These findings support a molecular link between monoaminergic antidepressant treatment and glutamatergic neurotransmission.

Cell type specific transcriptomic differences in depression show similar patterns between males and females but implicate distinct cell types and genes.

Major depressive disorder (MDD) is a common, heterogenous, and potentially serious psychiatric illness. Diverse brain cell types have been implicated in MDD etiology. Significant sexual differences exist in MDD clinical presentation and outcome, and recent evidence suggests different molecular bases for male and female MDD. We evaluated over 160,000 nuclei from 71 female and male donors, leveraging new and pre-existing single-nucleus RNA-sequencing data from the dorsolateral prefrontal cortex. Cell type specific transcriptome-wide threshold-free MDD-associated gene expression patterns were similar between the sexes, but significant differentially expressed genes (DEGs) diverged. Among 7 broad cell types and 41 clusters evaluated, microglia and parvalbumin interneurons contributed the most DEGs in females, while deep layer excitatory neurons, astrocytes, and oligodendrocyte precursors were the major contributors in males. Further, the Mic1 cluster with 38% of female DEGs and the ExN10_L46 cluster with 53% of male DEGs, stood out in the meta-analysis of both sexes.

Kabuki syndrome stem cell models reveal locus specificity of histone methyltransferase 2D (KMT2D/MLL4).

Kabuki syndrome is frequently caused by loss-of-function mutations in one allele of histone 3 lysine 4 (H3K4) methyltransferase KMT2D and is associated with problems in neurological, immunological and skeletal system development. We generated heterozygous KMT2D knockout and Kabuki patient-derived cell models to investigate the role of reduced dosage of KMT2D in stem cells. We discovered chromosomal locus-specific alterations in gene expression, specifically a 110 Kb region containing Synaptotagmin 3 (SYT3), C-Type Lectin Domain Containing 11A (CLEC11A), Chromosome 19 Open Reading Frame 81 (C19ORF81) and SH3 And Multiple Ankyrin Repeat Domains 1 (SHANK1), suggesting locus-specific targeting of KMT2D. Using whole genome histone methylation mapping, we confirmed locus-specific changes in H3K4 methylation patterning coincident with regional decreases in gene expression in Kabuki cell models. Significantly reduced H3K4 peaks aligned with regions of stem cell maps of H3K27 and H3K4 methylation suggesting KMT2D haploinsufficiency impact bivalent enhancers in stem cells. Preparing the genome for subsequent differentiation cues may be of significant importance for Kabuki-related genes. This work provides a new insight into the mechanism of action of an important gene in bone and brain development and may increase our understanding of a specific function of a human disease-relevant H3K4 methyltransferase family member.

Oxytocin receptor expression and epigenetic regulation in the anterior cingulate cortex of individuals with a history of severe childhood abuse.

Childhood abuse significantly increases the lifetime risk of negative mental health outcomes. The oxytocinergic system, which plays a role in complex social and emotional behaviors, has been shown to be sensitive to early-life experiences. While previous studies have investigated the relationship between early-life adversity and oxytocin, they did so with peripheral samples. We, therefore, aimed to characterize the relationship between early-life adversity and oxytocin receptor (OXTR) expression in the brain, using post-mortem human samples, as well as a rodent model of naturally occurring variation in early-life environment. Focusing on the dorsal anterior cingulate cortex, we compared OXTR expression and epigenetic regulation between MDD suicides with (N = 26) and without history of childhood abuse (N = 24), as well as psychiatrically healthy controls (N = 23). We also compared Oxtr expression in the cingulate cortex of adult rats raised by dams displaying high (N = 13) and low levels (N = 12) of licking and grooming (LG) behavior. Overall, our results indicate that childhood abuse associates with an upregulation of OXTR expression, and that similarly, this relationship is also observed in the cingulate cortex of adult rats raised by low-LG dams. Additionally, we found an effect of rs53576 genotype on expression, showing that carriers of the A variant also show upregulated OXTR expression. The effects of early-life adversity and rs53576 genotype on OXTR expression are, however, not explained by differences in DNA methylation within and around the MT region of the OXTR gene.

Association between the expression of lncRNA BASP-AS1 and volume of right hippocampal tail moderated by episode duration in major depressive disorder: a CAN-BIND 1 report.

The pathophysiology of major depressive disorder (MDD) encompasses an array of changes at molecular and neurobiological levels. As chronic stress promotes neurotoxicity there are alterations in the expression of genes and gene-regulatory molecules. The hippocampus is particularly sensitive to the effects of stress and its posterior volumes can deliver clinically valuable information about the outcomes of antidepressant treatment. In the present work, we analyzed individuals with MDD (N = 201) and healthy controls (HC = 104), as part of the CAN-BIND-1 study. We used magnetic resonance imaging (MRI) to measure hippocampal volumes, evaluated gene expression with RNA sequencing, and assessed DNA methylation with the (Infinium MethylationEpic Beadchip), in order to investigate the association between hippocampal volume and both RNA expression and DNA methylation. We identified 60 RNAs which were differentially expressed between groups. Of these, 21 displayed differential methylation, and seven displayed a correlation between methylation and expression. We found a negative association between expression of Brain Abundant Membrane Attached Signal Protein 1 antisense 1 RNA (BASP1-AS1) and right hippocampal tail volume in the MDD group (ß = -0.218, p = 0.021). There was a moderating effect of the duration of the current episode on the association between the expression of BASP1-AS1 and right hippocampal tail volume in the MDD group (ß = -0.48, 95% C.I. [-0.80, -0.16]. t = -2.95 p = 0.004). In conclusion, we found that overexpression of BASP1-AS1 was correlated with DNA methylation, and was negatively associated with right tail hippocampal volume in MDD.

Treatment-emergent and trajectory-based peripheral gene expression markers of antidepressant response.

Identifying biomarkers of antidepressant response may advance personalized treatment of major depressive disorder (MDD). We aimed to identify longitudinal changes in gene expression associated with response to antidepressants in a sample of MDD patients treated with escitalopram. Patients (N = 153) from the CAN-BIND-1 cohort were treated for 8 weeks, and depressive symptoms were assessed using the Montgomery-Åsberg Depression Rating Scale at 0, 2, 4, 6, and 8 weeks. We identified three groups of patients according to response status: early responders (22.9%), later responders (32.0%), and nonresponders (45.1%). RNA sequencing was performed in blood obtained at weeks 0, 2, and 8. RNA expression was modeled using growth models, and differences in the longitudinal changes in expression according to response were investigated using multiple regression models. The expression of RNAs related to response was investigated in the brains of depressed individuals, as well as in neuronal cells in vitro. We identified four RNAs (CERCAM, DARS-AS1, FAM228B, HBEGF) whose change over time was independently associated with a response status. For all except HBEGF, responders showed higher expression over time, compared to nonresponders. While the change in all RNAs differentiated early responders from nonresponders, changes in DARS-AS1 and HBEGF also differentiated later responders from nonresponders. Additionally, HBEGF was downregulated in the brains of depressed individuals, and increased in response to escitalopram treatment in vitro. In conclusion, using longitudinal assessments of gene expression, we provide insights into biological processes involved in the intermediate stages of escitalopram response, highlighting several genes with potential utility as biomarkers of antidepressant response.

Integrative DNA Methylation and Gene Expression Analysis in the Prefrontal Cortex of Mexicans Who Died by Suicide.

BACKGROUND: Suicide represents a major health concern, especially in developing countries. While many demographic risk factors have been proposed, the underlying molecular pathology of suicide remains poorly understood. A body of evidence suggests that aberrant DNA methylation and expression is involved. In this study, we examined DNA methylation profiles and concordant gene expression changes in the prefrontal cortex of Mexicans who died by suicide. METHODS: In collaboration with the coroner's office in Mexico City, brain samples of males who died by suicide (n = 35) and age-matched sudden death controls (n = 13) were collected. DNA and RNA were extracted from prefrontal cortex tissue and analyzed with the Infinium Methylation480k and the HumanHT-12 v4 Expression Beadchips, respectively. RESULTS: We report evidence of altered DNA methylation profiles at 4430 genomic regions together with 622 genes characterized by differential expression in cases vs controls. Seventy genes were found to have concordant methylation and expression changes. Metacore-enriched analysis identified 10 genes with biological relevance to psychiatric phenotypes and suicide (ADCY9, CRH, NFATC4, ABCC8, HMGA1, KAT2A, EPHA2, TRRAP, CD22, and CBLN1) and highlighted the association that ADCY9 has with various pathways, including signal transduction regulated by the cAMP-responsive element modulator, neurophysiological process regulated by the corticotrophin-releasing hormone, and synaptic plasticity. We therefore went on to validate the observed hypomethylation of ADCY9 in cases vs control through targeted bisulfite sequencing. CONCLUSION: Our study represents the first, to our knowledge, analysis of DNA methylation and gene expression associated with suicide in a Mexican population using postmortem brain, providing novel insights for convergent molecular alterations associated with suicide.

Non-CG methylation and multiple histone profiles associate child abuse with immune and small GTPase dysregulation.

Early-life adversity (ELA) is a major predictor of psychopathology, and is thought to increase lifetime risk by epigenetically regulating the genome. Here, focusing on the lateral amygdala, a major brain site for emotional homeostasis, we describe molecular cross-talk among multiple mechanisms of genomic regulation, including 6 histone marks and DNA methylation, and the transcriptome, in subjects with a history of ELA and controls. In the healthy brain tissue, we first uncover interactions between different histone marks and non-CG methylation in the CAC context. Additionally, we find that ELA associates with methylomic changes that are as frequent in the CAC as in the canonical CG context, while these two forms of plasticity occur in sharply distinct genomic regions, features, and chromatin states. Combining these multiple data indicates that immune-related and small GTPase signaling pathways are most consistently impaired in the amygdala of ELA individuals. Overall, this work provides insights into genomic brain regulation as a function of early-life experience.

Object Weighting: A New Clustering Approach to Deal with Outliers and Cluster Overlap in Computational Biology.

Considerable efforts have been made over the last decades to improve the robustness of clustering algorithms against noise features and outliers, known to be important sources of error in clustering. Outliers dominate the sum-of-the-squares calculations and generate cluster overlap, thus leading to unreliable clustering results. They can be particularly detrimental in computational biology, e.g., when determining the number of clusters in gene expression data related to cancer or when inferring phylogenetic trees and networks. While the issue of feature weighting has been studied in detail, no clustering methods using object weighting have been proposed yet. Here we describe a new general data partitioning method that includes an object-weighting step to assign higher weights to outliers and objects that cause cluster overlap. Different object weighting schemes, based on the Silhouette cluster validity index, the median and two intercluster distances, are defined. We compare our novel technique to a number of popular and efficient clustering algorithms, such as K-means, X-means, DAPC and Prediction Strength. In the presence of outliers and cluster overlap, our method largely outperforms X-means, DAPC and Prediction Strength as well as the K-means algorithm based on feature weighting.

Cocaine-related DNA methylation in caudate neurons alters 3D chromatin structure of the IRXA gene cluster.

Epigenetic mechanisms, like those involving DNA methylation, are thought to mediate the relationship between chronic cocaine dependence and molecular changes in addiction-related neurocircuitry, but have been understudied in human brain. We initially used reduced representation bisulfite sequencing (RRBS) to generate a methylome-wide profile of cocaine dependence in human post-mortem caudate tissue. We focused on the Iroquois Homeobox A (IRXA) gene cluster, where hypomethylation in exon 3 of IRX2 in neuronal nuclei was associated with cocaine dependence. We replicated this finding in an independent cohort and found similar results in the dorsal striatum from cocaine self-administering mice. Using epigenome editing and 3C assays, we demonstrated a causal relationship between methylation within the IRX2 gene body, CTCF protein binding, three-dimensional (3D) chromatin interaction, and gene expression. Together, these findings suggest that cocaine-related hypomethylation of IRX2 contributes to the development and maintenance of cocaine dependence through alterations in 3D chromatin structure in the caudate nucleus.

Stimulation of L-type calcium channels increases tyrosine hydroxylase and dopamine in ventral midbrain cells induced from somatic cells.

Making high-quality dopamine (DA)-producing cells for basic biological or small molecule screening studies is critical for the development of novel therapeutics for disorders of the ventral midbrain. Currently, many ventral midbrain assays have low signal-to-noise ratio due to low levels of cellular DA and the rate-limiting enzyme of DA synthesis, tyrosine hydroxylase (TH), hampering discovery efforts. Using intensively characterized ventral midbrain cells derived from human skin, which demonstrate calcium pacemaking activity and classical electrophysiological properties, we show that an L-type calcium agonist can significantly increase TH protein levels and DA content and release. Live calcium imaging suggests that it is the immediate influx of calcium occurring simultaneously in all cells that drives this effect. Genome-wide expression profiling suggests that L-type calcium channel stimulation has a significant effect on specific genes related to DA synthesis and affects expression of L-type calcium receptor subunits from the CACNA1 and CACNA2D families. Together, our findings provide an advance in the ability to increase DA and TH levels to improve the accuracy of disease modeling and small molecule screening for disorders of the ventral midbrain, including Parkinson's disease.

Integrated genome-wide methylation and expression analyses reveal functional predictors of response to antidepressants.

Major depressive disorder (MDD) is primarily treated with antidepressants, yet many patients fail to respond adequately, and identifying antidepressant response biomarkers is thus of clinical significance. Some hypothesis-driven investigations of epigenetic markers for treatment response have been previously made, but genome-wide approaches remain unexplored. Healthy participants (n = 112) and MDD patients (n = 211) between 18-60 years old were recruited for an 8-week trial of escitalopram treatment. Responders and non-responders were identified using differential Montgomery-Åsberg Depression Rating Scale scores before and after treatment. Genome-wide DNA methylation and gene expression analyses were assessed using the Infinium MethylationEPIC Beadchip and HumanHT-12 v4 Expression Beadchip, respectively, on pre-treatment peripheral blood DNA and RNA samples. Differentially methylated positions (DMPs) located in regions of differentially expressed genes between responders (n = 82) and non-responders (n = 95) were identified, and technically validated using a targeted sequencing approach. Three DMPs located in the genes CHN2 (cg23687322, p = 0.00043 and cg06926818, p = 0.0014) and JAK2 (cg08339825, p = 0.00021) were the most significantly associated with mRNA expression changes and subsequently validated. Replication was then conducted with non-responders (n = 76) and responders (n = 71) in an external cohort that underwent a similar antidepressant trial. One CHN2 site (cg06926818; p = 0.03) was successfully replicated. Our findings indicate that differential methylation at CpG sites upstream of the CHN2 and JAK2 TSS regions are possible peripheral predictors of antidepressant treatment response. Future studies can provide further insight on robustness of our candidate biomarkers, and greater characterization of functional components.

Mutations in ACTL6B Cause Neurodevelopmental Deficits and Epilepsy and Lead to Loss of Dendrites in Human Neurons.

We identified individuals with variations in ACTL6B, a component of the chromatin remodeling machinery including the BAF complex. Ten individuals harbored bi-allelic mutations and presented with global developmental delay, epileptic encephalopathy, and spasticity, and ten individuals with de novo heterozygous mutations displayed intellectual disability, ambulation deficits, severe language impairment, hypotonia, Rett-like stereotypies, and minor facial dysmorphisms (wide mouth, diastema, bulbous nose). Nine of these ten unrelated individuals had the identical de novo c.1027G>A (p.Gly343Arg) mutation. Human-derived neurons were generated that recaptured ACTL6B expression patterns in development from progenitor cell to post-mitotic neuron, validating the use of this model. Engineered knock-out of ACTL6B in wild-type human neurons resulted in profound deficits in dendrite development, a result recapitulated in two individuals with different bi-allelic mutations, and reversed on clonal genetic repair or exogenous expression of ACTL6B. Whole-transcriptome analyses and whole-genomic profiling of the BAF complex in wild-type and bi-allelic mutant ACTL6B neural progenitor cells and neurons revealed increased genomic binding of the BAF complex in ACTL6B mutants, with corresponding transcriptional changes in several genes including TPPP and FSCN1, suggesting that altered regulation of some cytoskeletal genes contribute to altered dendrite development. Assessment of bi-alleic and heterozygous ACTL6B mutations on an ACTL6B knock-out human background demonstrated that bi-allelic mutations mimic engineered deletion deficits while heterozygous mutations do not, suggesting that the former are loss of function and the latter are gain of function. These results reveal a role for ACTL6B in neurodevelopment and implicate another component of chromatin remodeling machinery in brain disease.